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We consider Bloch electrons on the honeycomb lattice under a uniform magnetic field with 2-Kp/q 
flux per cell. It is shown that the problem factorizes to two triangular lattices. Treating magnetic 
translations as Heisenberg-Weyl group and by the use of its irreducible representation on the space 
of theta functions, we find a nested set of Bethe equations, which determine the eigenstates and 
energy spectrum. The Bethe equations have simple form which allows to consider them further 
in the limit p, g — * oo by the technique of Thermodynamic Bethe Ansatz and analyze Hofstadter 
problem for the irrational flux. 
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I. INTRODUCTION 

The quantum Hall effect (QHE) on graphite thin films 
has been observed recentlyi^^^ . It has the honeycomb lat- 
tice with one tt electron per lattice site. The only Fermi 
levels {Ep = 0) are two points of the Brillouin zone where 
the conduction band and the valence bands touch forming 
cones. Low energy excitations have a linear spectrum as 
for the massless relativistic Dirac particles in 2+1 dimen- 
sional space-timei^i^. InfiiM the Hall effect on graphene 
on the basis of Dirac electrons was investigated. 

The problem of electrons in two dimensional periodic 
potential in a magnetic field has been attracted a lot of 
attention due to the unexpected nature of the commen- 
surability and frustrations.9d°,"d2.i3d4,i5d6ii7,i8,i9,2i,22^ 

It has let the existence of Cantor set energy spectra, and 
the wavefunctions strongly depend on wether the flux is 
commensurate or incommensurate with the lattice. It 
was observed^^^ that when the magnetic flux correspond- 
ing to unit cell of the periodic structure of the system 
is incommensurate with the potential, the one particle 
spectrum exhibits multifractal behavior like the Cantor 
set (see alsoi^). Then, after discovering the topological 
character of the conductance in the model by linking it 
with the Chern class in the linear bundle of the Bloch 
wave functioni^iM, its possible connection with quantum 
Hall effect become evident. 

Usually a solid with well localized atomic orbits is mod- 
elled by the lattice and the investigations were carried out 
by the use of a tight-binding Hamiltonian.. The pres- 
ence of zero modes was pointed out ini^ and the model 
with nearest-neighbor and next-to-nearest-neighbor hop- 
ping was studied inii. Wiegmann and Zabrodin^^ have 
observed the presence of a hidden quantum group in the 
Hofstadter problem on the regular lattice and initiated 
studies22i2ii2^ of the integrable structure in it. The ap- 
proach allowed to express the spectrum and the Bloch 
wave function at the mid-band points as solutions of the 
Bethe equations typical for completely integrable quan- 
tum systems. In the article^S the problem have been 
analyzed by the use of cyclic representations of quan- 



tum groups at any momenta, but the explicit form of 
the Bethe equations was not found(besides the midband 
point). 

IniSi, authors have used the fact that magnetic transla- 
tions in Hofstadter problem are forming Heisenberg-Weyl 
group and analyzed its irreducible representations on the 
basis of theta functions^''. This approach allowed us to 
find explicit form of the Bethe equations for the spectrum 
and the wave- functions not only at the midband point, 
but for any momenta. 



Motivated by the discovered unconventional Hall effect 
in graphene we study in this article the problem of Bloch 
electrons in a uniform magnetic field with flux 2i:p/q per 
cell on the honeycomb lattice. We show that the prob- 
lem possesses a chiral factorization and it reduces to two 
(i? and L) triangular sublattices. Following the line of^i, 
we look for the solution of the Schrodinger equation for 
energy and eigenfunctions in the space of irreducible rep- 
resentation of the Heisenberg-Weyl group, namely, in the 
space of theta functions with characteristics. By the use 
of nesting procedure, similar to the one appeared in Alge- 
braic Bethe Ansatz for some of integrable models (see for 
example^^iSi2&) , we have transformed Schrodinger equa- 
tion into the simple set of nested Bethe equations for 
the spectral parameters. We have found a spectrum of 
the model and eigenstates as a function of the solution 
of the nested Bethe Anzats equations, which depend on 
momenta k. 



The form of the nested Bethe equations allow to inves- 
tigate them in thermodynamic limit g ^ oo by the tech- 
nique developed by Takahashi^^ and Gaudin^^ for inte- 
grable models and called Thermodynamic Bethe Ansatz. 
This limit is important for the most interesting irrational 
flux case of the Hofstadter problem since q gives the num- 
ber of bands and we face a Cantor like behavior of the 
spectrum. Irrational number can be approached by the 
rationales p/q with p,q ^ oo. 
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II. ELECTRONS ON THE 2D HONEYCOMB 
LATTICE IN A UNIFORM MAGNETIC FILED 



Let us start with formulation of Hofstadter problem 
on the honeycomb lattice C (see FIG.l). The honeycomb 
lattice possesses chiral factorization, namely one can con- 
sider C — C'~ + as a sum of two triangular lattices £^ 
and each of which is isomorphic to the dual of the 
honeycomb lattice. Sites of the lattices and on 
the Fig^are marked with black and yellow dots respec- 
tively. The links of the honeycomb lattice are connecting 
the sites on with the three nearest neighbor sites on 




FIG. 2: Hopping vectors on the honeycomb (Ccr) and triangu- 
lar (/io-) lattices 
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FIG. 1: The representation of the honeycomb lattice as a joint 
of triangular (red) and (blue) lattices. 



The Hamiltonian is given by 



nG£«,o-=l,2,3 



d+ Cn), 



(1) 



where Cfj.d and dfjo-^* ,di,^ are the annihilation and 
creation operators of electrons at the sites n and n + 
of the triangular lattices and CJ^ respectively. 
^n,n+e„ — " ^n-(-e„.n IS the vcctor potential of a mag- 
netic field with the strength perpendicular to the plane 
of the lattice and ta^^ are hoping amplitudes between 
the nearest neighbors. We choose them as ti, t2 and 
according to directions of the vectors Co- (see Fig. 2). For 
the homogeneous model ti = t2 = ^3- 

In the following we consider the diagonalization prob- 
lem of in one-particle sector. In this sector the action 
of the translation operator by vectors ico- can be written 
as 



= Yl + (2) 



where the standard bra- and ket- vectors are used: 

\^) = ci|wac), {n\ = {vac\cn and \n + e^) = d't_^_gjvac) , 
{n + e„\ = {vac\dfi+g^. 

We consider the case when the magnetic flux per 
hexagon cell, 



exp(i$) = Y\_ exp{iAri,fri), 

hexagon 



(3) 



is rational: $ ~ 2tt^, where p and q are mutually prime 
integers. The product in 10) is performed on anticlock- 
wise direction. 

In terms of bra- and ket- vectors the Hamiltonian 
can be written as 



\fi){n + e'er I 



7=1,2,3 



+ is+e-„,se*'^"+'^-"|n + e^)(n|) +/I.C. (4) 

Let us now consider the square of the Hamiltonian Q). 
It is straightforward to sec from the expression that 
contain only hopping terms within triangular lattices 
and , namely 



where 



H 



R{L) 



E 



H'- = HB. + nL, (5) 



(6) 



with (see FIG. 2) 



^1 = 62 - 63, 

Ma 



M2 = ei - 63, 

ei - 6*2, 



(7) 



and 



t5 



n,n±pi — 6 ^ ^2^31 



in,n — ^1 + ^2 "I" ^3 



(8) 
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This decomposition reveal the hidden chiral structure 
of the Hofstadter problem on the honeycomb lattice re- 
ducing it to the problem on the triangular lattice. There- 
fore, in the remaining part of this article, we will consider 
only the triangular lattice. 

The problem of Bloch electrons on the triangular lat- 
tice and in uniform magnetic filed was considered recently 

:„29 



III. ELECTRONS ON THE TRIANGULAR 
LATTICE AND IN UNIFORM MAGNETIC FIELD 



Denote lattice sites by vectors n = nifli+n2fi2- Taking 
different gauges one can obtain various equivalent forms 
of (0). In the following we will use the Landau gauge, 
which is 



(9) 



Then the Hamiltonian © is invariant under translations 
S±il^^ S±p2 and Slp,^: 

So, the problem of diagonalization of Ti. reduced to its 
diagonalization on each eigenspace of S'^p^ , S'i^j and 

Sj-p^. The latter is a g-dimensional space ^'(A:), spanned 
by Bloch wave functions 



E 

ni = nmodq 



-ikn 



(10) 



which satisfy 

iln{k) = 1pn+q{k) 

S±p,Mk)=e^'''^Mk) 
Let us define also the generators of magnetic translations 

by 

T±^. =^e*^-'^±-|n)(n±A2a|, ^ = 1,2,3. (11) 

n 

It is easy to find out, that in Landau gauge © 

T±p,^e^^/'T^p,T±p,, (12) 

since the area of the triangles in honeycomb lattice 
(Fig.l) is the half of the area of the hexagon. 

The generators T±p-^ and T±p^ satisfy the following 
commutation relations 



and form Heisenberg-Weyl group. Here we used the no- 
tation 



, \ ( -P 

q = exp 1 *— ) = exp [ TTi — 



Note, that is any gauge Si - — Ti - (in the Landau 
gauge moreover 5*;;^ = T)rJ. 

The action of magnetic translations T±p-^ , T±p^ on 
Bloch functions 4'n(k) has the following form: 

T±p^Mk) = e^^'^'i^n^iik) 
T±pMk) = eT''=-±»*Vn(fc) (14) 

In case of translational invariant distribution of hopping 
parameters t±o- = in,n±p„, c = 1,2,3, and t = tfl^fj, the 
Hamiltonian can be written in terms of the generators 
of magnetic translations by 

iiaTp^+i-.T^P.) + T.^\n){n\- (15) 

cr=l,2,3 n 

Note that on '^{k) the g-th power of magnetic transla- 
tions are scalars: 



-^±Pl — ^ l*^' ^±p2 ~ ^ 



rpq 



(16) 



In order to solve the Schrodinger of equation TLip = Sip 
analiticaly we will follow the technique developed in^^ for 
the Hofstadter problem on the regular lattice. 



IV. HEISENBERG-WEYL GROUP. 
REDUCTION OF THE PROBLEM TO THE 
BETHE EQUATIONS 



Let us recall the representation of Heisenberg-Weyl 
group p3(l on the space of complex functions^^. It can 
be constructed in the following way. Define actions Sb 
and Ta, a,b E Z on the space of analytic functions on the 
complex plain as 



Sbf{z) = fiz + b), 
Ta{T)f{z) — exp{TTia'^T + 2'Kiaz)f{z + ar), 

for some t <E C. Then 



(17) 



T±p-^T-!^p^ q T-j-p^T-j-p-^ , 



(13) 



SaSb = Sa+b,Ta{T)%{T) = Ta+b{T) 

Sa%{T) = exp(2^ia6)rb(r)5a. (18) 

Consider the space of theta functions with characteris- 
tics e(«)(z,r) which are invariant with respect to the 
subgroup generated by S±q and 7±i(t). They form a q- 
dimensional space Thq(T) and have in the fundamental 
domain with vertices (0, T,q,q -\- r) precisely q zeroes. 
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The space Thq(r) forms an irreducible representation 
of Heisenberg-Weyl group H13|) generated by 

T±p, - a^^±e'*/2^±i(r)5±p. (19) 
One can choose the basis Qa.o{z, t) of Thg as follows: 
Qadz.r) = \%{t)Q\{z,t), 



where a = 0, . . . , - 6(z,t) = 9o(^,r) 

is the standard theta function and Oa_f,(z,T) = 
TaSbQ{z,T), a e modi, b e modi' {^) is the 

' N 



M 



notation for theta functions with characteristics 
In this basis 

'^iiOa.o = exp(±27rm)ea,o 



Comparing these equations with (|14|l and H16|l we obtain 
for parameters a±, I3± 

a± = exp(=Fifci) (3± = cxp(^jfc2), 



The space is identified with the space of theta func- 
tions with characteristic 



The space Thg(r) have some analog with the space 
of polynomials of degree q — 1. The polynomial de- 
composition analog for the 0'^*-' is the following. Every 
0(9) £ xhg can be represented in the form 



Z — Zi T 



e(«)(z,r)=7exp(-^)n0 



(20) 



where X]i=i Zi — kr + nq, k,n ^Z, < k,n < q, 
7 G C. Then for our Hamiltonian on triangular lattice 



+ e-**r_,,,)+tii2(e-'*r^3+e**r_,,3), (21) 



which was induced from the Hofstadter problem on the 
honeycomb lattice and the Schrodinger equation Tiip — 
£ip can be written in terms of theta functions. It is con- 
venient to write them in terms of 



„ , , , vriT ni . \ ^ / 1 — T 
ei(z,r) = exp ( — + y - TTiz 1 e I z H —,t 



1 / ' 



■!/'„(fc) ~ Q:i.fi{z,T) 



which obey 8i(0,t) = 0. After rescaling of the argu- 
ments of 0i(z, r): z —f qz, T ^ qr and performing addi- 
tional substitution: z — > z -I- we obtain 



-l)f/3+e^e-^nQi( 



i=l 



■ 2^ I - l + 2fc — q 



-(-l)^/3_e-t+^nQi(^ 

1=1 

n ei (z - z. - I , r) + «_e'f +-^+2... -Q 
i=i ^ 

^ne.(« 



z - Z.J H , r 
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t) = (£:- 3)^61(2 -z„T 



(22) 



where 



(23) 



r 



and £ is the eigenvalue of the Hamiltonian H21|) on the So, inserting the zeros z = z^, i = 1, . . . , g of Oi into the 
triangular lattice. left hand side we obtain the system of q equations for the 

The roots of right and left side of (|22|l must coincide. spectral parameters z^: 



5 



(-l)P/3+e**-'^ n ei (z, - + ^, r) + (-If /3_e-**+'^ [] Gi (z, - - ^, r) 



+ a+e 



i\Zj - Zi ,r +a_e 



l-2fc + ., _^2 



i^-\-7rtT — -—27rtZj'] 



1 



i=l 



Zj Zi -\- , T 



z, - z,; , r I = 



(24) 



This is an analog of Bethe Ansatz equation for the 
problem under consideration, but it has no convinient 
form for further investigations and needs to be simplified. 

By the use of the nesting procedure in the Appendix A 
we show that the Schrodinger equation H22|l transforms 
to the very simple form 



(£ - 3) - 7 



nLiei(z- 



.(12) 



73 



nLi0i(^ 



(25) 



where parameters 73 and 7 are defined by the formu- 
las (jASI - IA7II , while the set of Bethe equations (|^ is 
equivalent to the following larger, but much simpler set 
of nested Bethe equations 



nLi0i( 



.(1) 



nLi0i( 



.(1) 



,(2) 



q' I 



(26) 



nLiQiC 



nLi0i( 



.(2) 



,(3) 



a- 
a 



(27) 



I+p \ 



nLi0i( 



.(3) 



„i*(i+2fc)-27riTa^-47rjz<.-'> 



(28) 



,(12) 



.(1) 



= (-1)P+1 



nLi0i( . 
nLi0i(^i 



.(12) 



(2) N 
q ' ' 



^+nLi0i(^i 
nLi0i(^i 



.(2) 



.(1) 



nLi0i(^i 



2p 



(29) 



and 



-27ri2i ®l(^ 



J12) 



nLi0i(^. 
s+ nLl 01 (^1 



A3) 



r) 



^+nLi0i( 



Zl 



nLl 01 (^1 



(2) 



J12) 



J3) , (p+r) 



nLi0i(4 



(2) 



z«,r 



).ei(zi-z, + ^,r) 



(30) 



We have obtained a set of five equations for the set of five 



l,2,---g. The pa- 



parameters Zi, zl"'"'' , Zj-^"* , zp'' , zl^^^ , 
rameters a±, f3± and 5± depend on momenta k. There- 
fore the solutions Zi, zf'\ zf'\ zf\ zf''^\ i — l,2,---q of 
the equations l|26|) - (|3UII depend on the momenta too. 

Both terms in the sum on the left hand side of H25() 
are elliptic functions, i.e. 2-periodic, in the fundamental 
domain with vertexes (0, r, 1, 1 -I- r), and each of them 
has precisely q zeroes and q poles there. Let us suggest 
that all poles are single. It is known that such an elliptic 
function can be written in the form 



g 



/(z)=C+> A,;C(Z-Z0, 



1=1 



(31) 



where ^(z) is a ^-function of Weierstrass for the same 
periodicity domain. Formula (|31|l is an analog of decom- 
position of rational functions with single poles in terms 
of functions 

Z — Z j 

Applying (|31|l like decomposition for elliptic functions 
to the expression (|25() for the energy one can verify that 
the conditions for the annulation of the residues of the 
poles precisely are coinciding with H26|l - (|30|l (or equiva- 
lently H24|) '). Therefore the energy E is independent of z 
constant function C in this decomposition and one can 
find (see formulas (|A9|) and ljA6|l in the Appendix) 

£ = 3 + (_i)P5^e**(3-fe)+^»^^^^^+2'^«^i (32) 
nLl ei(zi - + r) ■ 61 (zj^^) - zf ), r) 

'nLl 01(^1-4'^ + ^ 



,r) • 61 (z(^')-z„r) 



This expression defines the spectrum of the states 
in the Hofstadter problem on the honeycomb lat- 
tice as a function of momenta because the solutions 
Zi, zf'\ zf'\ zf \ zl^'^\ j = 1, 2, ••• of the equations 
(|26I30|) are momentum dependent. 

This procedure reminds Nested Bethe Ansatz where 
we have a nested chain of Bethe equations (see as an 
example the set of nested Bethe equations in t-J— and 
staggered t-J modelsS^). 

Although we have increased the amount of Bethe equa- 
tions (|26|I - H3U|) together with the amount of parameters 
to be defined, each of the equations (|26H3UI) now is much 
simpler and allows to apply the technique of Thermo- 
dynamic Bethe AnsataSLSS for their investigation. After 
taking logarithm of the left and right hand sides of these 
equations, one can reduce them to the set of integral 
equations in the most interesting limit of q oo, nec- 
essary in the irrational flux case. This fact justifies the 
procedure of nested Bcthc Ansatz. 

Nested Bethe equations ^ EO)) now have a form, 
which allows to investigate them further by use of Ther- 
modynamic Bethe Ansatz. One can take logarithm of 
the left and right hand sides and transform it to integral 
equation in the thermodynamic limit g — > oo. 
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APPENDIX A: NESTING OF BETHE ANSATZ 
EQUATIONS 



Since any term in the linear space of functions O^"?) E 
Thq can be represented as a product of theta functions 
(|20|1 with some particular choice of zeros at Zj , j — 
l,2,---q, the sum of the first and second products of 
theta functions in the expression (|22|l can be represented 
as 



■ * I 27rifcp 

+ (-lf/3_e-*6+— 



ne 



i=l 



p 

1 

Q 



l\Z - Zi - -,T 



7in®i 

i=l 



(1) 



(Al) 



1=1 



= 72 



nei(z-zf,r), (A2) 



(-l)P(5+e*3+"^— 
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■ JjBlfz - Zi + -,T 



i=l 

29 , _,_l + 2fc-g 



q 



-27Tiz+- 



l\Z- Zi 



T + p 



,(3) 



,r) (A3) 



1=1 



Right hand sides of the equations IjAip , HA2|I and (|A3|) 

u 4- (1) (2) (3) • -, o 

become zero at z = , z — z^ , z — Zj , j — 1, 2 ■ ■ ■ q, 
respectively, which imposes restrictions on the parame- 
ters Zi, z'l^\ z^-^\ zf^^i — 1, 2, • • • g in the form of Bethe 
equations 

The parameters 71 , 72 and 73 are z independent con- 
stants, as follows from the arguments presented in the 
previous section by use of representation of elliptic func- 
tions via Weierstrass C, function (|31|l . Therefore they can 
be fixed by inserting any convenient value z into the cor- 
responding equation. By putting values z = zi+p/q^ z = 
zi — r/g, z = zi+pI q+r/ q into the equations (jAl|) . (|A2p 
and 1A3I) respectively, which annulate one of terms in the 
left hand sides, one can solve this equations for 7's and 
obtain 



(A4) 



71 


= (-1; 




-fe) 




nt 


.i0i(^i - 








i0i(^i- 




72 


= a-|_e 




q ^ 






.161(^1- 






nt 


i0i(^i- 




73 


= (-1; 




-k)-\-7riT- — -\-2Trizi 




nt 


.i0i(^i - 






nL 


,ei(zi- 





(A5) 



(A6) 



Consequently, the sum of the third and forth terms in 
(respectively the sum of fifth and sixth terms) define 



ne 



1 ( z - Zi - - , r 



Now let us continue further this procedure, which re- 
minds nesting Bethe Ansatz in the some of integrable 
model9^^i2SiS&, and transform the sum of right hand sides 
of (jAljl and (|A2|) into the product of theta functions with 
the zeros at the new positions z'f''^^ 
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£ to the very simple form 



(2) N 



1=1 



1=1 



(A7) 



(2) 

Now again, one can insert the value z = z\ into the 
equation above and obtain 



7 71 



nLi0i(4 



nLi0i(-f'-4"\r)' 



(A8) 



where 71 defined by the expression in (|A4|) . The spectral 
f 12) 

parameters are the zeros of the left(and, therefore, 
of the right) hand side of the equation (jA7|) . Hence, by 
use of (jA4p and (|A8p we obtain Bcthe equations H29|) for 

With this reformulations of the sums of theta functions 
we bring Schrodinger equation H22fl for the eigenenergy 



{£~i)\{Q,(z-z,,T)=-t\{Q,{z~z^^\T) 



+l:,X{Qi{z-zf\T) (A9) 



i=l 



('12') 

By inserting z — z\ into this expression (the first 
term on the right hand side becomes zero) and using 
(|A6(I for 73 one can obtain simple expression l|32|) for 
the energy. This complets Bethe Ansatz solution of the 
Hofstadter problem on the honeycomb lattice. 

The set of nested Bethe equations H26l - I3()|l . are equiv- 
alent to the single set (|24ll . We have more equations in 
(|26l - l30ll than in H24|l . but their advantage is in their sim- 
plicity. Each of mentioned equations contains only two 
product of theta functions, which will allow to implement 
the technique of Thermodynamic Bethe Ansata2L2& for 
further investigations. 
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